Quantum Computing in Arrays 
Coupled by 'Always On' Interactions 



S. C. Benjamin 1,2 and S. Bose 3 

1 Ctr. for Quantum Computation, Clarendon Lab., Univ. of Oxford, OXl 3PU, UK. 
2 Dept. of Materials, Parks Road, Univ. of Oxford, 0X1 3PH, UK. 
3 Dept. of Physics and Astronomy, University College London, 
Gower St., London WC1E 6BT, UK. 

It has recently been shown that one can perform quantum computation in a 
Heisenberg chain in which the interactions are 'always on', provided that one 
can abruptly tune the Zeeman energies of the individual (pseudo-) spins. Here 
we provide a more complete analysis of this scheme, including several general- 
izations. We generalize the interaction to an anisotropic form (incorporating the 
O ■ XY, or Forster, interaction as a limit), providing a proof that a chain coupled 
in this fashion tends to an effective Ising chain in the limit of far off-resonant 
spins. We derive the primitive two-qubit gate that results from exploiting abrupt 



o 
o 

(N 
03 



a: 



Zeeman tuning with such an interaction. We also demonstrate, via numerical 
simulation, that the same basic scheme functions in the case of smoothly shifted 
Zeeman energies. We conclude with some remarks regarding generalisations to 



two- and three-dimensional arrays. 



There has recently been considerable interest in the question of whether one can per- 
form quantum computation (QC) in Heisenberg-type systems (e.g. interacting electron 
spins) when the interaction is 'always-on'^ 3 -. This question follows on from a work con- 
cerning Heisenberg systems in which the interactions are presumed to be switchable, either 
individually 4,5 6 or collectively^. Numerous proposals exis t 8 ' 9 ' 10 for experimental realization 
of such a model, however interaction switching is liable to prove very challenging to real- 
ize, and this motivates the interest in 'always-on' interactions. In Ref. [l] we proposed a 
scheme for exploiting a simple one-dimensional Heisenberg chain with constant, isotropic 
nearest neighbor interactions. The scheme involved adjusting the single-spin level splittings 
(the Zeeman energies) to bring neighbors in and out of resonance with one another. We 
exploited the fact that far off-resonance spins do not exchange energy, but rather interact 
in an Ising 'ZZ' form. We argued that by separating the qubit-bearing spins by passive 
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'barrier' spins, one can negate this residual interaction (thus achieving a passive state for 
the array) - yet one can invoke an interaction on demand simply by bringing a barrier into 
resonance with its neighbors. 

In the present paper we elaborate on several aspects of that earlier Letter, and we provide 
certain extensions. Whereas previously we considered only one specific form for the inter- 
action, i.e. the isotropic Heisenberg form (a x cr x + a Y o Y + a z a z ), we now generalize our 
arguments to accommodate different magnitudes for the in-plane and perpendicular compo- 
nents. Thus we subsume the prior isotropic form, and the purely planar "XY" interaction, 
as special cases. There is a wide variety of promising physical systems associated with this 
family of interactions (for the isotropic limit, see e.g. Refs. |s IqI1iq|] . and for the anisotropic 



case, Refs. 
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13]). The XY limit is also referred to as the Forster interaction, espe- 



cially when studied in the context of excitonic exchange in biological molecules. With this 
generalized form of interaction, we first present an analysis of the effect of far off-resonant 
neighbors in a long chain, obtaining the anticipated Ising-like form as the lowest order term. 
We then explain in detail how Zeeman tuning can be exploited perform an elementary two- 
qubit gate, and we show how the resulting unitary operation depends on the Z versus XY 
asymmetry in the interaction. 

Whereas the original paper assumed a perfectly abrupt transition between on-resonant 
and far off-resonant Zeeman energies, here we follow our analysis with a numerical simulation 
demonstrating that smoothly changing Zeeman energies can implement the gate process 
equally well. This observation considerably increases the practicality of the scheme. Finally, 
we discuss the generalization to two- and three-dimensional arrays. 



Analysis of Heisenberg Chain with Large Zeeman Discrepancies 

The analysis is presented in full in Appendix I. Here we summarize it. We start from a 
total Hamiltonian H given by 

H H single ~t~ -"intj 

where 

H single ^ ] ^j^i " 
3 

and the exchange interaction is as follows, where the factor a allows for a possible anisotropy 
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between the in-plane and z-direction components. 

H in t = J^ a 3 a 3+l + °i a 3+l + a(X 3 a 3+l) 
3 

= i(E Vj°7+i + a i a t+i) + Ja E °f°f+i> 

where a" 11 = <r x ± ia y . Here and below, the sum ranges over all N qubits, but subscripts 
such as j + 1 are understood to be modulo N, i.e. we assume a closed circular topology. 
This considerably simplifies the analysis, but it is not a real constraint - in the limit of large 
chains the open and closed topologies will be equivalent. 
We rewrite H = Hi + H 2 where 

H 1 = 52B j <rf + Ja'E<TftTf +1 . 

3 3 

and 

3 

Notice that H\ is simply the Hamiltonian for an Ising spin chain with varying Zeeman 
energies. We will find that this term dominates the time evolution when the spins are far 
off-resonance with their neighbors; the contribution of H 2 then vanishes. 

Our approach is to exploit the Trotter formula to manipulate the time evolution operator 
into a form that can be recognized as Ising and non-Ising parts. This is detailed in Appendix 
I. The exact expression for the time evolution is found to be: 

U(t) = R(t) exp (-mt) where R = {lC=i exp =fH R (f) 

Hr(v) = E X 3^)° j + vj+i + X}(ri)ajaf +1 with ^(77) = exp(i7 7 (A j + 2Ja(af +2 - af_ ± )) 
3 

where Aj = 2(B j+ i — Bj). The right hand term in U(t) is the pure Ising chain evolution 
we seek, but the subsequent 'residual' operator R is more complex. In the second part of 
Appendix II we expand R as a power series and inspect the terms. We conclude that, for a 
regular chain with a characteristic A (such as an ABABAB.. chain where Aj = (— 1) J A), 
the time evolution can be written as 

U(t) = (1 - 8P(t)) exp (-Wit) 

where 5 = J/A, for some finite operator P(t) whose magnitude does not increase with A. 
Thus for any given time period t the non-Ising evolution will be negligible if A 
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is sufficiently large compared to J. Assuming that we can dynamically change a Aj, 
switching it between zero and a large value, we can then exploit this result to produce a 
form of 'gate' for quantum computation. 



Exploitation of the Heisenberg-to-Ising Transition to Perform QC 

Assume that we have some array in which every pair of adjacent spins is far off resonance 
from one another, i.e. Aj ^> J, Vj. Now assume that we abruptly tune one (or more) of the 
spin Zeeman energies so that we have a triplet ABA where energies A and B are comparable. 
Let us refer to these spins by the labels 1 to 3, and similarly label the external neighboring 
spins as and 4. Suppose spins 0, 2 and 4 are initially in state | \). Since spin remains 
far off resonance from 1, their interaction is effectively of the Ising form Jaa^af . Similarly 
the interaction between 3 and 4 is Jaa^af. Moreover, those external spins (having only 
an Ising interaction with their neighbors) are 'frozen' in the | |) state thus their interaction 
with the triplet reduces to Jacrf and Jacrf , and the dynamics of the triplet are described 
by the Hamiltonian: 

-^triplet — H 7 „ 



'zeeman -^int 
#zeeman = (A + aJ)(<rf + (if ) + £<jf 



#int = J E °f a f+ 1 + 



J 3 a j+l 



aafaf +1 



3=1,2 



In the following we will use the notation J X y = J, Jz = a J, a = A + J z (the effective 
Zeeman energy of spins 1 and 3) and b = B for consistency. The Hamiltonian is easy to 
analyze; the states | fit) and | J, J, I) of course remain eigenstates while the remaining states 
form two distinct subspaces. For the 'up' subspace spanned by {| |Tt)> I TIT)) I TTI)} we 
have Hamiltonian and eigenvectors given by 



Hu = 61 + 2 J XY 
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With corresponding energies = b, E± = b + Jxy(p ± S p ). Here p = (a — b — J z )/ Jxy 



and S p = y/8 + p 2 . Similarly for the complimentary 'down' space {| Til), I ITI)> I I IT)} we 
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have 



H D = -bI + 2J XY 
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With energies E% = -b, E± = b+J X Y(p±S p ), where q = (b-a-J z )/JxY and S q = + q 1 ■ 
Now, we know that the initial computational qubit states are 



|00> = 


lit!) 


<= 


composed of 


\+) D and \-) D 


|01> = 


1 ITT) 


<= 


composed of 


\a)u, \+)u and I — > tr 


|10> = 


ITU) 


<= 


composed of 


\a)u, \+)u and \-) v 


lll> = 


ITTT) 


<= 


eigenstate 





During the gate operation, the states (other than |11)) will rotate within their subspaces. 
We must arrange to 'revive' both the 1 00) state and the states |01) & 1 10} at the same 
instant, i.e. we must arrange that at some time tR the central spin is in the definite state 
| |) for all computational basis states. (Note that this condition does permit a net rotation 
in the plane defined by 1 01} & 1 10} ) . Thus at that moment we can effectively switch off the 
exchange interaction (by switching to far off-resonant Zeeman energies) and we will have 
performed some unitary transform in the computational basis. Whether such a transform 
constitutes a useful gate depends on entanglement criteria as mentioned later. The times 
for which 1 00) revives are determined by — E®. The times at which a state, initially in 
the 1 01) , 1 10) plane, returns to that plane are determined by E+ — E*i. Now the parameter 
which we can experimentally vary is the Zeeman detuning a — b; although there may be 
various detunings for which the revivals coincide (which could be found numerically), there 
is one value that is immediately obvious by inspection: a — b = (corresponding to tuning 
the central barrier spin to A + J z) ■ In this case we see that p — q — — Jz/ Jxyi 5p — Sq — 
^8 + (Jz/ Jxy) 2 and thus both revivals coincide at time tR = irh(8Jx Y + J§)~^ . At this 
instant, the transformation in the computational basis (|00), 1 01) , 1 10) , |11)} is given by the 
following matrix (neglecting a global phase) 

/n n n n \ 



u = 



10 

iQs Qc 
Qc iQs 
W 



Here Q = — exp(i<f>), s/c = sin / cos(0) and W = — exp(— 2i<p) with = ~ (SJ^y/Jf + 
The phases in this matrix are with respect to the passive state of the device (i.e. if we 
had not tuned the triplet into resonance), under the assumption that the resonance was 
achieved by shifting the Zeeman energy of the central spin. (If in fact the Zeeman energies 
of the qubit-bearing spins were adjusted to achieve resonance, then we simply have the 
above matrix together with two trivial single qubit Z gates.) This transformation U is 
entangling, and is therefore adequate to construct a universal gate set when combined with 
single qubit gates^. Using the procedure described in Refs. |14l)l5j | one can confirm that 
no more than four uses of this gate are required to form a Control-NOT, for a wide range 
of Jz including the Jz = and J = Jz cases, which represent the XY interaction and the 
isotropic Heisenberg interaction, respectively. It is easier to appreciate the nature of the 
transform if we apply a couple of single-qubit Z-rotations; defining 

/ 

Z{0) = 



exp(z#) 
exp(— it 



then neglecting a global phase, 



Z x ty).Z 2 {$)-U 
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-Q's iQ'c 

iQ'c -Q's 

1 



(1) 



Uj 



\ 



(2) 



Here ip = |(1 — (8J 2 /J| + 1)~^) and Q' = Q 2 while s/c are as before. Notice that for the 
Jz = limit, i.e. the case of a pure XY interaction, then the primitive matrix U takes a 
particularly simple forrn^ 

/ i 
0-10 
0-10 
-1 

using which one can construct a CNOT with only two applications, as shown in Fig. 1(b). 
In this limit, the dressed matrix (0) is recognizable as the "iSWAP" which has been studied 
in the context of an XY interaction between adjacent qubits 17 . Indeed, in the limit of a 
strict XY interaction, one might choose to abandon the barrier spin architecture completely, 
and adopt a trivial architecture in which qubits are adjacent (since the primary function of 
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R v_=^(^z>-CJ Y ) R z =^(I+ia Y ) 

FIG. 1: (a) Schematic showing the basic two-qubit gate. Letters A, B, € denote Zeeman energies; 
A and B are fixed but € is abruptly changed as shown by the graph to the right. The operation is 
completed when the central spin 'revives' into the state | |). The resulting primitive two-qubit gate 
is entangling for all values of the anisotropy parameter a = Jz/Jxy- In the special case of Jz = 
the gate has a particularly simple form (eqn. |2J) and a corresponding circuit for the control-NOT 
operation (b) requires only two such primitives. 

the barrier spins is to negate the effect of the residual Ising interaction, absent for the pure 
XY form). 

Note the second form of gate presented in Ref.flJ can also be generalised to anisotropic 
Heisenberg interactions, although it does require some finite Z component since this is 
exploited to accumulate a phase during the gate operation. 

In Ref. |l( and in the above analysis, we consider an abrupt change from far off-resonance 
spins into resonance. This may be difficult to achieve in many otherwise promising imple- 
mentations, therefore it we now investigate the effect of smooth switching. Figure 2 shows 
various profiles for the dynamically changing Zeeman energy of the central spin, given that 
the Zeeman energies of the outer spins are static (Fig. 2(a) corresponds to analytic treatment 
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FIG. 2: (a) An abrupt Zeeman shift, corresponding to our analytic treatment, (b) & (c) Numerical 
simulation demonstrating that other, smooth functions can also suffice. In (b) we use a 'gentle' 
switching function of the form cos 2 ( 7r ^j^ L ) > where to is the time at which the switching process 
begins and t& = 1.25 is the switch duration. After a 'flat' period, the function is reversed to return 
us to the large 'passive' detuning. In (c) we construct a 4 th power sinusoidal function for a sharper 
profile 18 . In each case we find that simply choosing the correct duration for the central constant 
phase allows complete^ barrier revival. In these graphs, the time axis is in units of %/ Jxy, an d 
we arbitrarily chose an anisotropic interaction with Jz/Jxy = 0.7. 

given above). For numerical convenience we have built these switching profiles as piece-wise 
combinations of analytic functions, as denned in the figure caption. In both cases (b) and 
(c) we fixed the time for the switching transition to the arbitrary choice £a = 1-25 and varied 
just a single parameter, the time for which the detuning is zero. The values shown in the 
Figure provided a complete revival of the central spin for all qubit basis states, just as in 
the case of the abrupt transition. The specific transformation achieved in the qubit basis 
(i.e. the analogue of eqn. ([TJ) is of course different for these smooth switching profiles, but 
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FIG. 3: Structures that are efficient in terms of Rq, the ratio of number of qubits stored to total 
number of spins. In one dimension, using the gate process analysed here (c.f. Fig. 1) the value of 
Rq is necessarily 1/2. For higher dimensions, if one employs the same gate process (by 'inserting' 
a unique barrier between each pair of adjacent qubits) then the best value of Rq for a regular 
structure is 2/5; structures meeting this limit are shown in (b) and (c). If, however, one allows 
each barrier to separate multiple qubits, then higher ratios are possible - for (d) Rq = 3/5. 

it remains strongly entangling and therefore equally suitable as a primitive two-qubit gate. 

The analysis presented in the present paper has been phrased in terms of a one- 
dimensional array (a). However, the basic gate construction, involving two qubit bearing 
spins and one barrier spin, can immediately be generalised to many geometries in either 
two, or three dimensions. In principle one can produce a suitable structure by taking any 
arrangement of qubit-bearing spins, and introducing a barrier spin between each (hitherto) 
adjacent pair. One possible measure of the efficiency of the implementation would be the 
ratio of qubit-bearing spins to total number of spins, which we can denote Rq. The value 
Rq = | corresponds to the one-dimensional arrangement (Fig 3(a)). For for a two, or 
higher, dimensional geometry at least some of the qubits must of course have three or more 
neighbors. If we restrict ourselves to considering regular structures in which every qubit has 
the same number of neighbors, then it is apparent that the highest possible value of Rq is 
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2/5. Two arrangements which achieve this value are the hexagonal geometry Fig. 3(b), and 
the 3D structure illustrated in Fig. 3(c). 

In order to do better than this ratio it would be necessary for barrier spins to do 'double 
duty' in the sense that each barrier could not be unique to a specific qubit pair. Figure (d) 
shows an example arrangement achieving Rq = 3/5 in 2D. Note (d) is the compliment of 
(b), i.e. the qubit and barrier roles are reversed; similarly, one could reverse Fig. 3(c) for 
a 3D form. In such a structure, bringing a barrier into resonance with its neigbors would 
initiate a three-qubit gate process - to successfully complete the gate one would require the 
simultaneous revival of all qubit basis states at some subsequent moment. As the number 
of qubits involved increases, this quickly becomes infeasible (see Appendix II), but both the 
three qubit gate shown in Fig. 3(d), and a four qubit variant, do appear possible 21 . Of 
course, such multi-qubit gates are quite exotic and may be rather inefficient primitives for 
implementing algorithms. 

In conclusion, we have extended the results presented in Ref. in several significant 
respects. The first, fundamental generalization is from a pure isotropic Heisenberg interac- 
tion to a more general anisotropic interaction, including the in-plane "XY" interaction as a 
special case. All the results presented here incorporate this generality. We have provided a 
proof that an interaction of this general form tends to a simple Ising interaction in the limit 
of far off-resonance neighbors. We have presented an analysis of the basic gate of Ref. [l] 
in with this general interaction, and exhibited the resulting primitive two-qubit gate. In 
the special case of an XY interaction, we note that the gate has an especially simple form 
and we provided an explicit circuit for an efficient CNOT based on this primitive. We also 
consider the effect of a non-abrupt switching of the Zeeman energy: by numerical simula- 
tion we demonstrate that simply varying the duration of the on-resonance phase (while the 
switching time remains constant) allows one to achieve the necessary revival of the barrier 
spins, and therefore abrupt switching is not a requirement of the scheme. Finally we have 
remarked upon the simplicity of generalizing to two- and three-dimensional arrays, noting 
that the array geometry then determines the scheme's cost in terms of the proportion of 
barrier spins. 

SCB wishes to acknowledge support from a Royal Society URF, and from the Foresight 
LINK project "Nanoelectronics at the Quantum Edge". 
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Appendix I: Analysis of Heisenberg Chain with Large Zeeman Discrepancies 



Given the definition of H = Hi + H2 introduced in the main body of the paper, we can 
proceed to use the Trotter formula to write the time evolution operator U(t) as 

U(t) — ^exp (— iHit/n) exp (— iH 2 t/n) n as n — > 00. (3) 

Now we will seek to move all Hi terms to the right, thus separating the Ising and non-Ising 
parts. Note first that since 

Wi , <Tj I = Wi ><Tj<Tk\ =Q \. a i °3 ' a k ° m \ = 

we can write the following, using r = t/n, 

exp(-iHir)= JJexp(-i.B,-T<7f) Y[(exp(-iaJraf af +1 ) (4) 
3 3 

and in fact we can reorder these terms as we wish. Moreover we can use 

exp(irBjaf) = cos(rBj)l + i sm(rBj)af 
exp(irJaaf af +1 ) = cos(r J a)l + ism(r J a) a? cr? +1 (5) 

We will also find it useful to employ 

~z„± 1+ + z _,_ + . „z + „±„Z ta\ 

where a ± are as defined in the main body of the paper. We will introduce a generalisation 
of H 2 , 

Hf = ^(E WjataT+i + Wjarati). 

where the Wj are any functions involving scalar constants and <rf for any/all k. Now expand 

00 1 —i T 

eM-iH^r) = Ezj,^-(E Wrt*r +1 + Wjaja^r * (7) 
and note the following using (j3J) and © 

expi-irBjcrf ) [ W j <Tf<rJ +1 + W}oja+ +1} 

= { W 3 a t a 7+i + W k7 a t+i) exp{2irB j af)exp{-irB j af) 
= ^e-^Wja+aJ^ + e 2irB 'Wjaj <r+ +1 ) exp(-iS i rof ) 
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and similarly 

exp(-trB J+1 af +1 )(W j( x+aJ +1 + W}aja+ +l ) 

= (e 2iTB ^W j( J+a]- +1 + e-^+WV:^) exp(-iB j+1 raf). 

Then 

3 

= (exp{irA k )W k a^a k+1 + exp(-iTA k )W j k a k a^ +1 ) TJ(exp (-irBjcrf ) 

where Aj = 2{Bj + \ — Bj). Now combining this with (@| and ((7J) we can write 
exp(— iH\t) exp(-iH^r) 

= H(exp(-taJrafaf +1 ) exp(-iH^r) J] ex^-iBjraf ) (8) 

where Vj = exp(irAj)Wj. Now we can commute the remaining left side product through to 
the right in a similar way. Again using (jSJ) and © we note that: 

= yWj<Jj~<j] +1 + Wjaj'<Tj +1 exp^iraJaf^af) 

= , W j a t a 7+^ + W l a 7 a t+i\ ex P(2?™Jo-^_ 1 o- J z )exp(-?raJcr^_ 1 (j J z ) 



= ^jexp{—2iTaJa z _ l )a^Oj +1 + WjexY>{2iraJoj_ l )o : j Cj+i. exp(— zraJcr^cr 
Similarly 

exp{-iraJa? +1 af +2 ) ^a+aj^ + W}aja+ +l} 

= ^Wj exp(2iraJaf +2 )a^aJ +1 + Wjexp(—2iraJaf +2 )a^aj +1 ^ exp(—iraJaf +1 af +2 ) 

However, for the afaf +1 term we see that 

[exp(-iraJafaf +1 ) , Wa^aj +l + W^ajaf +1 ] = 

since there is a double sign inversion. Then combining these three results we can write 

J[exp(-iraJafaf +1 ) ) { W k a+a k+1 + W 1 k a k a+ +1 ) 

= W k exp(2ir Ja(af +2 - <T%_d)<T$ 'a k+1 
+W ] k exp(-2irJa(a^ +2 - of^))^ <r£ +1 T[exp(-iraJaf af +1 ) 
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Now combining this with (JJJ) and (JBJ) we have 

exp(— iH\r) exp(— iH^r) = exp(— iH^r) exp(— iH\r) 

where Qj = Wjexp(ir(Aj + 2Ja(crf +2 - crf_i)). 

Now because these Qj fit within the original definition of Wj (i.e. they are simply "func- 
tions involving scalar constants and erf for any/all fc"), we can just repeat the argument to 
commute all terms exp(zfZir) to the far left. The term originally identified as the m th ele- 
ment exp(iH2r) in the Trotter expansion (J3J) will have m terms "exp(— if/ir)" pass 'through' 
it, and will thus accumulate a final Q(m) = exp(im^(Afc + 2Ja(a^ +2 — cjf-i))- So the exact 
expression for the time evolution finally becomes: 

U(t) = R(t) exp {-iHxt) where R = {UZ=i e*P } ri ^ 00 

#*(»7) = E X i(^>7+i + X i^) (T 7 (7 /+i with = exp(^(A, + 2Ja(af +2 - af_,)) 

j 

The right hand term in U(t) is the pure Ising chain evolution we seek, but the subsequent 
'residual' operator R is more complex. We would like to show that it tends to unity as 
5j = J/Aj — > for all j. Now we cannot simply integrate the terms in the product R 
since they do not commute, and thus we cannot immediately gather the elements with a - 
coefficient. Therefore we proceed by making the expansion exp A = 1 + A + A 2 /2 + ...: 

if exp {-iH R {— )t/n) = if 1 - iH R (—)t/n + l/2(iH R (— )t/n)f + ... 

as n — > oo. We cannot truncate this series since is not small, but we will seek to gather 
and sum all terms of given order in H R . We introduce 

n t 

9a = 22 -H R {mb/n) 

m=a " 

f drj H R {rf) 

J rj=at/n 

we can evaluate indefinite integrals 20 as follows. Defining A = A , Pj = A/Aj and 5 = J/A, 
^ / ^i(^) = -^-a^-fa) ~ / drjXj(rj) = l —XjX}{rj) 
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where 

_ - 8g»g(l + af_ 1 crf +2 ))(l - H«2 ~ *f-i)) 
Xj " 1 - 16a 2 5j ~ Pj 

with the approximation holding in the limit that all 5j = ^- <C 1. Note that and thus 

Xj, is a modest ratio in our periodic chains (e.g. for an ABAB.. chain pj = (— i)- 5 ; for an 

ABC ABC. chain might run 1, 1, — |, 1, 1, — |, .. say). We can write the indefinite integral 

= jK(r,) (9) 

Using KQ defined above we can write 

g a = ^(K(t)-K(at/n)) (10) 

Now returning to the expansion, the lowest order in H R is of course 1, and the sum of all 
terms of 1 st order in H R is precisely —igo = \5{K{t) — K(0))}. Thus so far we are seeing 
the anticipated behavior: the 'residual' part of the dynamics, after the Ising-like behavior 
is allowed for, appears to vanish with 5. However, since we are using an expansion in H R (), 
where H R () is not small, we should evaluate and sum the higher terms. Let us use the 
symbol Sn to represent the sum of terms of order H R Q ; then we have already found 
Si = \8{K{t) - K(0)}, and 

t n 1 

S 2 = (-i-) 2 J2{-(H R (mt/n)) 2 +H R ( m t/n)J2Mpt/2)} 

n m =l ^ p>m 
t U 1 

= -R J E{- (W»)) J + W n ) E H R (pt/2)} 

77 — Z — 

m=l p>m 

now the factor {t/n) 2 causes the first term here to vanish in the limit 77 — > 00, since it 
contains only 77 terms each of Order(H R () ~ J). For the second term 

S 3 = - f H R ( V ) f H R (rj)d V d V (11) 
J Ut| ) 

but the inner integral is given by (jl(J|) so that 

S2 = ~Y Io Hr{v) t K{t) ~ K{r]) ) dr] (12) 

This integral can be fully evaluated^ but the key point is that it can already be seen to be of 
order 5 (or less). Note that the 5 = ^ factor cannot be absorbed by the remaining integral 
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since the variable A occurs only as a phase ~ exp(iAr). Thus in expanding and evaluating 
the integral we will see some terms with an additional factor of ^, and in the special case 
that a term exhibits cancellation of the A elements in the phase we would apply a factor of 
order unity - but we can never introduce a factor of A. 

Generalizing this observation we can consider Sn- This involves terms of the form 
H R (mit/n)H R (m 2 t/n)...H R (rriNt/n) for some set of integers m\ > m 2 > ...m N . By the 
same reasoning above, we can neglect terms where two or more of the m ; are the same 
value, since they collectively constitute a negligible portion ^ of the sum as n — > oo. Then 
we find 

Sn = {~i) N f f H R {i 2 ) f .... f H R {i N ) d& d£ 2 ...d£ N 

JO J$l J§2 J£N-1 

(j) 1 *- 1 ft ft r-t r-t 

= { -^—6 / Hnfa) / H R (£ 2 ) / .... / (K(t) - Kfo-i)) dfc d& ...d^N-i 

2 Jo J£l J& J^n-2 

And as before we can argue that although the remaining N—l integrals may produce 
additional factors of 1/A, they cannot absorb any. Thus the factor 5 will remain and we can 
conclude that all terms in the expansion S n (n > 1) are of order 5 or less. Therefore 
the time evolution operator is 

U(t) = (1 - 5P(t)) exp (-iHit) 

for some finite operator P(t) whose magnitude does not increase with A. This is the result 
presented in the main body of the paper. 

Appendix II: Regarding Revivals 

In the discussion of two and three dimensional arrays, we stated that it will be difficult 
to achieve the crucial simultaneous 'revivals' for multi-qubit gates involving more than a 
few qubits. Of course, one can observe that if we choose any detuning A — B for which 
the revival periods of the various qubit basis states are related by irrational factors (i.e. 
the general case), then there will eventually be a complete revival to any desired accuracy 
(although never perfect). However one would typically need to wait an extremely long time 
for the level of precision required for QC and therefore this type of revival is not a practical 
choice. Instead we seek to arrange rapid revivals by looking for values of the detuning (and 
potentially, other parameters) such that the various revival periods are all related by small 
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rational factors. Fulfilling this condition will become unfeasible as the number of qubits 
increases. 
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